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Abstract 

We improve the theoretical estimates of the critical exponents for the three- 
dimensional Heisenberg universality class. We find 7 = 1.3960(9), v = 
0.7112(5), T] = 0.0375(5), a = -0.1336(15), (3 = 0.3689(3), and 5 = 4.783(3). 
We consider an improved lattice c/)^ Hamiltonian with suppressed leading scal- 
ing corrections. Our results are obtained by combining Monte Carlo simula- 
tions based on finite-size scaling methods and high-temperature expansions. 
The critical exponents are computed from high-temperature expansions spe- 
cialized to the (j)^ improved model. By the same technique we determine 
the coefficients of the small-magnetization expansion of the equation of state. 
This expansion is extended analytically by means of approximate parametric 
representations, obtaining the equation of state in the whole critical region. 
We also determine a number of universal amplitude ratios. 

PACS Numbers: 75.10.Hk, 75.10.-b, 05.70.Jk, 11.15.Me 
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TABLE I. Recent experimental estimates of the critical exponents for Heisenberg systems. 





Material 


7 


/5 


6 


Ref. |] (1980) 


Ni 




0.354(14) 




Ref. § (1981) 


Fe 




0.367(5) 




Ref. I] (1995) 


Ni 


1.345(10) 


0.395(10) 


4.35(6) 


Ref. § (1995) 


GdaBrC 


1.392(8) 


0.365(5) 


4.80(25) 


Ref. § (1995) 


GdalC 


1.370(8) 


0.375(8) 


4.68(25) 


Ref. § (1999) 


Tl2Mn207 


1.31(5) 


0.44(6) 


4.65(15) 


Ref. § (2000) 


Lao.82Cao.i8Mn03 




0.383(9) 




Ref. § (2000) 


Lao.gsCao.osMnOs 


1.39(5) 


0.36(7) 


4.75(15) 


Ref. |] (2000) 


Gd(OOOl) 




0.376(15) 




Ref. llOl (2000) 


Gd2Cu04 


1.32(2) 


0.34(1) 




Ref. Ill (2000) 


CsoPdao (liq) 


1.42(5) 






Ref. lul (2000) 


CsoPdao (sol) 


1.40(8) 






Ref. |1| (2001) 


GdS 




0.38(2) 




Ref. PI (2001) 


CrOa 


1.43(1) 


0.371(5) 




Ref. f^ (2001) 


Lao.gCao.aMnOs 


1.45 


0.36 





I. INTRODUCTION AND SUMMARY 

According to the universality hypothesis, some features of continuous phase transitions — 
for instance, critical exponents and scaling functions — do not depend on the microscopic 
details of the systems, but only on few global properties, such as the space dimensionality, 
the range of interaction, and the symmetry of the order parameter. These features define a 
universality class. In this paper, we consider the three-dimensional Heisenberg universality 
class, which is characterized by a three-component order parameter, 0(3) symmetry, and 
short-range interactions. 

The Heisenberg universality class describes [0] the critical behavior of isotropic magnets, 
for instance the Curie transition in isotropic ferromagnets such as Ni and EuO, and of 
antiferromagnets such as RbMnFa at the Neel transition point. In Table | we report some 
recent experimental results. It is not a complete review of published results, but it is useful 
to get an overview of the experimental state of the art. In the table we have also included 
results for the well-studied doped manganese perovskites Lai_j.A2;Mn03, although the nature 



of the ferromagnetic transition in these compounds is still unclear [|I5|. The Heisenberg 
universality class also describes isotropic magnets with quenched disorder. Indeed, since 
a < 0, the Harris criterion states that disorder is an irrelevant perturbation. The only effect 
is to introduce a correction-to-scaling term jtl^'^'^ with Ajis = —a. The experimental results 
confirm the theoretical analysis [|T6|, as it can be seen from Table y (older experimental 
results with a critical discussion are reported in Ref. P3). The prediction for Adis has been 



checked in perturbative field theory |^ and experimentally p6|JT8|JT9[] . 

Beside the exponents 7, /3, and 6 there are also a few estimates of the specific-heat 
exponent a, in most of the cases obtained from resistivity measurements: a ~ —0.10 in Fe 



and Ni |2^; a = -0.12(2) in EuO [|28|; a = -0.11(1) in Fe^Nigo-^BigSi 13; a = -0.11(1) 
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TABLE II. Recent experimental estimates of the critical exponents for Heisenberg systems 
with quenched disorder. 







Material 


7 


/3 


S 


Ref. 11?! (1994) 


FeioNiroBiigSi 


1.387(12) 


0.378(15) 


4.50(5) 


Ref. 1 


17| (1994) 


FeisNieyBiigSi 


1.386(12) 


0.367(15) 


4.50(5) 


Ref. 1 


17| (1994) 


Fei6Ni64Bii9Si 


1.386(14) 


0.360(15) 


4.86(4) 


Refs. 


|T8|,|T9ll (1995) 

.1 II II V / 


Fe2oNi6oPi4B6 


1.386(10) 


0.367(10) 


4.77(5) 


Refs. 


p|,ITl (1995) 


Fe4oNi4oPi4B6 


1.385(10) 


0.364(5) 


4.79(5) 


Ref. 1 


2| (1997) 


FegiZrg 


1.383(4) 


0.366(4) 


4.75(5) 


Ref. 


20| (1997) 


FesgCoZrio 


1.385(5) 


0.368(6) 


4.80(4) 


Ref. 


2| (1997) 


Fe88Co2Zrio 


1.389(6) 


0.363(5) 


4.81(5) 


Ref. 


2| (1997) 


Fe84Co6Zrio 


1.386(6) 


0.370(5) 


4.84(5) 


Ref. 


2T| (1999) 


Fe1.85Mn1.15Si 


1.543(20) 


0.408(60) 


4.74(7) 


Ref. 


2l| (1999) 


Fei.5oMni.5oSi 


1.274(60) 


0.383(10) 


4.45(19) 


Ref. 


2| (2000) 


Fe86Mn4Zrio 


1.381 


0.361 




Ref. 


22| (2000) 


Fe82Mn8Zrio 


1.367 


0.363 




Ref. 1 


23| (2001) 


Fe84Mn6Zrio 


1.37(3) 


0.359 


4.81(4) 


Ref. 1 


23| (2001) 


Fe74Mni6Zrio 


1.39(5) 


0.361 


4.86(3) 



in RbMnFg [|| . 

Aim of this paper is to substantially improve the precision of the theoretical estimates 
of the critical exponents. For this purpose, we consider an improved lattice Hamiltonian 
that is characterized by the fact that the leading correction to scaling is (approximately) 
absent in the expansion of any observable near the critical point. Moreover, we combine 
Monte Carlo (MC) simulations and analyses of high-temperature (HT) series. We exploit 
the effectiveness of MC simulations and finite-size scaling (FSS) techniques to determine 
the critical temperature and the parameters of the improved Hamiltonians [^0|-|37|, and 
the effectiveness of HT methods to determine the critical exponents for improved models, 
especially when a precise estimate of the critical point is available. This approach has 
already been applied to the three-dimensional Ising and XY [0,^] universality classes, 
achieving a substantial improvement of the estimates of the universal quantities that describe 
the critical behavior, such as the critical exponents and the scaling equation of state. 

We consider a simple cubic lattice and the nearest-neighbor (p^ lattice Hamiltonian 



(1) 



where (p^ is a three-component field. As shown in Ref. the Hamiltonian (|I]) is improved 
for A = A* ~ 4.4(7). Here, we extend the simulations of Ref. [^, obtaining a more accurate 
estimate of A*, A* = 4.6(4), and precise estimates of the critical /3c for several values of 
A. The analysis of the MC FSS results obtained for the improved 0^ lattice Hamiltonian 



already provides precise estimates of the critical exponents. As shown in Refs. p8|,|39|,|36 



an additional increase in precision can be obtained by combining improved Hamiltonians 
and HT methods. For this purpose, by using the linked-cluster expansion technique, we 
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computed HT expansions of several quantities and analyzed them using the MC results for 
A* and f3c- The final results significantly improve those obtained from the MC simulation. 
Moreover, they substantially improve those obtained using longer (21 orders) series for the 
standard Heisenberg model ||4^ . 

In Table we report our results for the critical exponents. We give the estimates 
obtained from the analysis of the MC data alone and those obtained by combining MC and 
HT techniques — they are denoted by MC+IHT, where the "I" refers to the fact that we 
are considering an improved model. The exponent a can be derived using the hyperscaling 
relation a = 2 — Su, obtaining a = —0.1336(15). We would like to stress that the good 
agreement between the MC and HT estimates is not trivial, since the critical exponents 
are determined from different quantities and limits. Indeed, the MC estimates are obtained 
from the analysis of the finite-size behavior for the size L — cx) at the critical point f3 = f3c, 
while the HT results are derived from the singular behavior of infinite-volume quantities as 

In Table |T| we compare our results with the most precise theoretical estimates that have 
been obtained in recent years. A more complete list of results can be found in Ref. |^6[. The 
results we quote have been obtained by Monte Carlo simulations (MC), from the analysis 
of the HT series for the standard Heisenberg model (HT), or by field-theory methods (FT). 
The MC results were obtained by applying FSS techniques to different Hamiltonians. Refs. 

model, 
The HT re- 



4^-§4| studied the standard 0(3)- vector Heisenberg model, Ref. |33] the improved 
and Ref. an isotropic ferromagnet with double-exchange interactions 



suits of Ref. were obtained analyzing the 21st-order HT expansions for the standard 
0(3)-vector model on the simple cubic (sc) and on the body-centered cubic (bcc) lattice. 
The FT results of Refs. |]46|-51| were derived by analyzing perturbative expansions in dif- 
ferent frameworks: fixed- dimension expansion (6th- and 7th-order series, see Refs. ) , 
e-expansion (to O(e^), see Refs. p9| , |60| ), and {d — 2)-expansion (to 0[{d — 2)'^], see Refs. 
pT|^53[). We quote two errors for the results of Ref. [|^: the first one (in parentheses) is 
the resummation error, and the second one (in brackets) takes into account the uncertainty 
of the fixed-point value of the coupling, which was estimated to be approximately 1% 
in Ref. |^8|. To estimate the second error we use the results of Ref. where the depen- 



dence of the exponents on g* is given. The results of Ref. |^ were obtained by using the 
so-called scaling-field method (SFM). Refs. [p^ - |55| , |6l| present results obtained by approx- 
imately solving continuous renormalization-group (CRG) equations for the average action, 
which is approximated to lowest and first order of the derivative expansion. We also men- 
tion the HT results of Ref. |6^: they performed a direct determination of the exponent a 
obtaining a = —0.11(2), —0.13(2) on the sc and bcc lattice. Ref. [^| computes the critical 
exponents for a Heisenberg fluid by a canonical-ensemble simulation. Depending on the anal- 
ysis method, they find l/u = 1.40(1), 1.31(1), (3/u = 0.54(2), 0.52(1), and 7/1/ = 1.90(3), 
1.87(3). Overall, all estimates are in substantial agreement with our MC+IHT results. We 
only note the apparent discrepancies with the MC estimates of r] of Refs. ||4^j43[] , and with 
the FT results of Ref. However, the reliability of the error bars reported in Ref. |5S 



is unclear: indeed, Ref. |^ analyzes the same perturbative series and reports much more 
cautious error estimates. 

We also present a detailed study of the equation of state. We first consider its expansion 
in terms of the magnetization in the HT phase. The coefficients of this expansion are directly 
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TABLE III. Estimates of the critical exponents. See text for the explanation of the symbols in 
the second column. We indicate with an asterisk (*) the estimates that have been obtained using 
the relations 7 = (2 — t])^, 2/3 = + ??), 6{1 + t]) = 5 — r]. 



Ref. 


Method 


7 


V 




/5 


b 


this work 


MC+IHT 


1.3960(9) 


0.7112(5) 


0.0375(5) 


0.3689(3)* 


4.783(3)* 


this work 


MC 


1.3957(22)* 


0.7113(11) 


0.0378(6) 


0.3691(6)* 


4.781(3)* 


H (2000) 


MC 


1.393(4)* 


0.710(2) 


0.0380(10) 


0.3685(11)* 


4.780(6)* 


H (2000) 


MC 


1.3909(30) 


0.6949(38) 




0.3535(30) 




@ (1996) 


MC 


1.396(3)* 


0.7128(14) 


0.0413(16) 


0.3711(9)* 


4.762(9)* 


11 (1993) 


MC 


1.389(14)* 


0.704(6) 


0.027(2) 


0.362(3)* 


4.842(11)* 


@ (1991) 


MC 


1.390(23)* 


0.706(9) 


0.031(7) 


0.364(5)* 


4.82(4)* 


@ (1997) 


HT sc 


1.406(3) 


0.716(2) 


0.036(7)* 


0.3710(13)* 


4.79(4)* 


(1997) 


HT bcc 


1.402(3) 


0.714(2) 


0.036(7)* 


0.3700(13)* 


4.79(4)* 


11 (1993) 


HT 


1.40(1) 


0.712(10) 


0.03(3)* 


0.368(6)* 




S (2001) 


FT d = 3 exp 


1.3882(10) 


0.7062(7) 


0.0350(8) 


0.3655(5)* 


4.797(5)* 


(1998) 


FT d = 3 exp 


1.3895(50) 


0.7073(35) 


0.0355(25) 


0.3662(25) 


4.794(14) 


II (1991) 


FT d = 3 exp 


1.3926(13) [39] 


0.7096(8) [22] 


0.0374(4) 






© (1977) 


FT d = 3 exp 


1.386(4) 


0.705(3) 


0.033(4) 


0.3645(25) 


4.808(22) 


(1998) 


FT e-exp 


1.382(9) 


0.7045(55) 


0.0375(45) 


0.3655(35) 


4.783(25) 


(1998) 


FT e-exp 


1.39* 


0.708 


0.037 


0.367* 


4.786* 


(2000) 


FT {d - 2)-exp 




0.695(10) 








(1984) 


SFM 


1.40(3) 


0.715(20) 


0.044(7) 


0.373(11) 


4.75(4)* 


(2001) 


CRG 




0.74 


0.038 


0.37 


4.78 


(2001) 


CRG 


1.374 


0.704 


0.049 


0.369 


4.720 


(1996) 


CRG 


1.465 


0.747 


0.038 


0.388 


4.78 



related to the zero-momentum n-point renormalized couplings, which are determined by an- 
alyzing their HT expansion. These results are used to construct parametric representations 
of the critical equation of state which are valid in the whole critical region, satisfy the correct 
analytic properties (Griffiths' analyticity), and take into account the Goldstone singularities 
at the coexistence curve. From our approximate representations of the equation of state we 
derive estimates of several universal amplitude ratios. Moreover, we present several results 
and different forms of the equation of state that can be compared directly with experiments. 
In particular, we can compare with the experimental results of Refs. p!7|j20| J^ , finding good 
agreement. 

The paper is organized as follows. In Sec. y we present our MC results. In Sec. |IT1| we 
present our results for the critical exponents obtained from the analysis of the HT series for 
the improved Hamiltonian (0). The equation of state is discussed in Sec. 0. We determine 
the small-magnetization expansion coefficients in Sec. |IV A| , give an approximate parametric 
representation of the equation of state in Sees. [IV B| and |1V C| , compute several amplitude 
ratios in Sec. [[VD| , and compare the theoretical results with experimental data in Sec. [IV E| . 
Details are reported in the Appendices. In App. ^ we present the analysis of the MC results 
and in App. |B| the analysis of the HT series. The expressions of several amplitude ratios in 
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TABLE IV. Final results for f3c and R* from fits with ansatz (^). In parentheses we give the 
statistical error and in brackets the error due to the corrections to scaling. 



R 


Za/ Zp 




t/4 




R* 


0.1944(1)[4] 


0.5644(1) [2] 


1.1394(1) [2] 


1.4202(2) [10] 




0.6862390(10) [12] 


0.6862386(11) [6] 


0.6862365(17) [12] 


0.6862369(17) [19] 



terms of the parametric representations are reported in App. y. 



II. MONTE CARLO SIMULATIONS 

The present MC simulations extend those of Ref. |p^. Here, we have considerably 



enlarged the statistics and added larger lattice sizes. Moreover, we have considered an 
additional quantity in order to improve the control over systematic errors. This way, we can 
increase the accuracy of A* and give precise estimates of the critical Pc for three values of A 
in a neighborhood of A*. For a detailed discussion of our methods, see Ref. [^. Details are 
reported in App. 0. 

We simulated the 0(3)-symmetric (p'^ model (|lD at A = 4.0, 4.5, and 5.0 on a simple cu- 
bic lattice with linear extension L in all directions. We measured the Binder parameter f/4, 
its sixth-order generalization f/g, the second-moment correlation length ^2nd, and the ratio 
Za/Zp, where Za is the partition function with anti-periodic boundary conditions in one of 
the three directions and Zp the corresponding one with periodic boundary conditions in all di- 
rections. The number of iterations for each lattice size and value of A was approximately 10^ 
for L = 6, 7, 8, 9, 10, 11, 12, 14, 16, 18, 20, 22, approximately 10^ for L = 24, 28, 32, 36, 40, 48, 
and 1-4x10^ for L = 56,64,80,96. With respect to Ref. [^, we have added new lattice 



sizes for all three values of A and considerably increased the statistics. In total, the whole 
study took about four years on a single 450 MHz Pentium HI CPU. 

In the first step of the analysis, we compute f3c and the fixed-point value of the dimen- 
sionless ratios R* for A = 4.5, using the standard cumulant crossing method of Binder. In 
particular, we fit our data with the ansatz 

R* = R{L,f3,), (2) 

where R* and Pc are free parameters. Our results are reported in Table |V[ Note that the 
four results for Pc are consistent within error bars. The statistical error of (3c obtained from 
Za/Zp and ^2nd/L is considerably smaller than that from U4 and Uq. As our final estimate 
we take Pc = 0.6862385(20), which is consistent with all four results. 

In addition, we determine jSc for A = 4.0 and A = 5.0. For this purpose, we use the ansatz 
(0), fixing L = 96 and taking the values of R* from Table |V[ Our results are summarized 
in Table 0. For both values of A, the results obtained from the four different choices of R* 
are consistent within error bars. As our final result we take that obtained from Za/Zp, since 
it has the smallest statistical error. 

Then, we locate A* by studying the scaling corrections to a quantity R defined in terms 
of two dimensionless ratios -Ri and R2. To define R, we fix a number Ri j which should be a 
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TABLE V. Results for /3c at A = 4.0 and 5.0 using only L = 96 and the ansatz R{/3c) = R* 



where R* is taken from Table IV. In parentheses we give the statistical error and in brackets the 



error due to the uncertainty on R*. 



X 


Za/ Zp 


S,2nd/L 


f/4 




4.0 


0.6843895(20) [15] 


0.6843887(21) [14] 


0.6843898(31) [20] 


0.6843898(31) [26] 


5.0 


0.6875638(21) [16] 


0.6875633(26) [15] 


0.6875655(34) [20] 


0.6875646(34) [26] 



good approximation to Rl, see Ref. |^6|]. Then, for a given value of A and L, we determine 
Pf{L, A) from 

RimX,(3f) = Rij . (3) 

In our analysis, /?/ is determined by taking either {Za/Zp)f = 0.1944 or (C2nd/-^)/ = 0.5644. 
Note that Pf approaches Pc as 

f3f = p, + CjL-'/'' + ..., (4) 

where the prefactor Cf depends on the choice of Rij. In particular, if Rij = R^, then 
Cf = and the leading corrections are proportional to L^'^/'^-'^, 
Next, we define R by 

R{L,X) = R2{L,X,f3f) . (5) 

Here, we take either or Uq as R2- Below, we often refer to R as i?2 at Rij- Up to 
subleading corrections, R behaves as 

R{L,X) ^R* + c{X)L-'^ . (6) 

The optimal value A* is obtained by solving c(A) = 0. We obtain A* = 4.6(4), 4.7(8), 4.7(8) 
and 4.6(8) from at {Za/Zp)f = 0.1944, at (6nd/^)/ = 0.5644, f/g at (^a/^p)/ = 0.1944 
and Uq at (^2nd/-^)/ = 0.5644, respectively. As our final result we quote 

A* = 4.6(4) (7) 

from at (^a/^p)/ = 0.1944. 

Finally, we compute the critical exponents u and t] using standard FSS methods. Usually, 
the exponent u is computed from the slope of a dimensionless ratio R at f3c- Here, following 
Ref. [|^, we replace (3c by which simplifies the error analysis, and determine u from the 
relation 



dR 
'8(3 



= a L^l\ (8) 

/3/ 



We study the derivative of all four quantities U^, f/g, Gnd/-^, and Za/Zp, and fix (3f by using 
either (^2nd/-^)/ = 0.5644 or {Za/Zp)f = 0.1944. We arrive at the final estimate 

i/ = 0.7113(11), (9) 
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where the error includes both the statistical and the systematic uncertainty. 

The exponent t] is computed from the finite-size behavior of the magnetic susceptibility: 

X\,, = cL'-^ . (10) 
In addition, we also use a fit ansatz that includes a constant background term: 

x\^^ = cL'-'^ + b. (11) 

As before, we fix /3/ by setting either (^2nd/-^)/ = 0.5644 or {Za/Zp)f = 0.1944. Our final 
MC estimate is 

rj = 0.0378(6). (12) 



III. CRITICAL EXPONENTS FROM THE IMPROVED HIGH-TEMPERATURE 

EXPANSION 



As shown in the case of the Ising and XY universality classes ||39| , |36| , the analysis 
of HT expansions for improved Hamiltonians with suppressed leading scaling corrections 
leads to considerably precise results even for moderately long series. In the present paper, 
the analysis of 20th-order HT expansions for the improved 0"^ lattice Hamiltonian, i.e. for 
A ^ A* = 4.6(4), allows us to substantially improve the accuracy of the estimates of the 
critical exponents. As we shall see, the results turn out to be more precise than those 
obtained in the preceding Section. They also significantly improve those obtained from the 
analysis of longer series (21 orders) for the standard Heisenberg model (which is recovered 



in the limit A — > oo) on the cubic and bcc lattices In this Section we report the results 
of our analyses of the HT series. The details are reported in App. 0. 

We determine 7 and u from the analysis of the HT expansion to 0(/9^°) of the magnetic 
susceptibility and of the second-moment correlation length. In App. p 2| we report some 
details and intermediate results so that the reader can judge the quality of our results 
without the need of performing his own analysis. They should give an idea of the reliability 
of our estimates and of the meaning of the errors we quote, which depend on many somewhat 
arbitrary choices and are therefore partially subjective. 

We analyze the HT series by means of integral approximants (lA's) of first, second, 
and third order. The most precise results are obtained biasing the value of f3c with its 
MC estimate. We consider several sets of biased lA's and for each of them we obtain 
estimates of the critical exponents. These results are reported in App. [B 2| . All sets of lA's 
give substantially consistent results. Moreover, the results are also stable with respect to the 
number of terms of the series, so that there is no need to perform problematic extrapolations 
in the number of terms in order to obtain the final estimates. The error due to the uncertainty 
on A* is estimated by considering the variation of the results when changing the values of A. 

Using the results reported in App. p 2| for the analysis at A = 4.0, 4.5, and 5.0, we obtain 



7 = 1.39582(10)[18] + 0.0015(A - 4.5), (13) 
u = 0.71111(5)[8] + 0.0009(A - 4.5). (14) 
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The number between parentheses is basically the spread of the approximants at A = 4.5 
using the central value of /3c, while the number between brackets gives the systematic error 
due to the uncertainty on Pc- Eqs. (|13|) and ( p!^ show also the dependence of the results on 
the chosen value of A. The A-dependence is estimated by using the results for A = 4.0 and 
A = 5.0. 

Using the MC estimate A* = 4.6(4), we obtain 

7 = 1.39597(10)[18]{60}, (15) 
i/ = 0.71120(5) [8]{36}, (16) 

where the error due to the uncertainty on A* is reported between braces. Thus, our final 
estimates are 

7 = 1.3960(9), (17) 

z/ = 0.7112(5), (18) 

where the uncertainty is estimated by summing the three errors reported above. 

Using the above-reported results for 7 and i' and the scaling relation 7 = (2 — 77)1/, we 
obtain t] = 0.037(2), where the error is estimated by considering the errors on 7 and u as 
independent, which is of course not true. We can obtain an estimate of r] with a smaller, 
yet reliable, error by applying the so-called critical-point renormalization method |^ to the 
series of x and This method provides an estimate for the combination rju. Proceeding 
as before, we obtain 

rju = 0.02665(18) + 0.00035(A - 4.5). (19) 

Taking into account that A* = 4.6(4), we find 

r/i/ = 0.02669(18) [14], (20) 

where the first error is related to the spread of the lA's and the second one to the uncertainty 
on A*, evaluated as before. Thus, 

= 0.0375(3) [2]. (21) 

Moreover, using the scaling relations, one obtains 

a = 2 -3z/ = -0.1336(15), (22) 
6 = = 4.783(3), (23) 



l + r] 

V 

2 

where the error of (5 has been estimated by considering the errors of v and 77 as independent. 



/3=^ (1 + 77) = 0.3689(3), (24) 
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IV. THE CRITICAL EQUATION OF STATE 



In this section we determine the critical equation of state characterizing the Heisenberg 
universahty class. The critical equation of state relates the thermodynamical quantities in 
the neighborhood of the critical temperature, in both phases. It is usually written in the 
form 

H ={B,)-'MM'-'f{x), (25) 
x = t{M/B)-'^/^, (26) 

where f{x) is a universal scaling function normalized in such a way that /(—I) = and 
/(O) = 1, and Be and B are the amplitudes of the magnetization on the critical isotherm 
and on the coexistence curve, 

M = B^H^I^ t = 0, (27) 

M = B{-tf H = 0, t<0. (28) 

Griffiths' analyticity implies that f{x) is regular everywhere for x > — 1. It has a regular 
expansion in powers of x, 

oo 

f{x) = l + J2fnX^ (29) 
and a large-x expansion of the form 

f{x)=X^J2fn^-'''^- (30) 

Moreover, at the coexistence curve, i.e. for x ^ —1 pS|-[72[ 

f{x)^Cf{l + xy. (31) 

The nature of the corrections to the leading behavior at the coexistence curve is less clear, 
see, e.g., Refs. [[70|- [73| , |56[| . From the scaling function /(x) one may derive many interesting 
universal amplitude ratios involving zero-momentum quantities, such as specific heat, mag- 
netic susceptibility, etc.... For example, the universal ratio Uq of the specific-heat amplitudes 
in the two phases can be written as (see, e.g., Ref. jT^ ) 



n=l 



oo 



^'=A=- ^ ^^^^ 
where, in the Heisenberg case for which — 1 < a < 0, 

= "^'"^'""'{^"^ + _ Ixr^fix) + r dy \yr' [f'{y) - /'(O) - y/"(0)] . (33) 

a - 1 a Jo 

We mention that the critical equation of state for the A^-vector model has been computed 
to O(e^) in the framework of the e-expansion and to 0{1/N) in the framework of the 



1/A^ expansion |p8 
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As our starting point for the determination of the critical equation of state, we determine 
the first few nontrivial coefficients of its small-magnetization expansion, by analyzing the 
corresponding HT series for the improved 0^ Hamiltonian. These results are then used to 
construct parametric representations of the critical equation of state which are valid in the 
whole critical region. Then, from our approximate representations of the equation of state 
we derive estimates of several universal amplitude ratios. This method has been already 
applied to the Ising universality class in three | |75|J5B[ | and two dimensions |]7^, and to the 



three-dimensional XY universality class |178| , |36 |. 



A. Small- magnetization expansion of the Helmholtz free energy 



We write the Helmholtz free energy as [76, 7£ 



Tfl 

= ^(M) - ^reg(M) = — A{z), (34) 

where m = 1/^, ^ is the second- moment correlation length, is the zero-momentum four- 
point coupling, and 

z = k\M\t-^ (35) 

where k is an appropriate amplitude ratio. The small-magnetization expansion of the free 
energy corresponds to the small-z expansion of A{z), 

^w=r'+r'+2:(^-^.-'^ (36) 

which also fixes the normalization k of z. Correspondingly, we obtain for the equation of 
state 

dM \M\ ' '• ^ ' 

with F{z) = dA/dz. Because of Griffiths' analyticity, J-'{M) has a regular expansion in 
powers oit for \M\ fixed. Therefore, F{z) has the large-z expansion 

F(z) =/^Ffc°°;2-'=/^. (38) 
fe=0 

The function F{z) is defined only for t > 0. For t < the equation of state is expressed 
in terms of a different function which is however analytically related to F{z) since the free 
energy and the equation of state are analytic on the critical isotherm t = for H 0. The 
two functions f{x) and F{z) are clearly related: 

z-'F{z) = F^fix), z = zox-^, (39) 

where Zq = kB. 
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TABLE VI. Estimates of (74, rg, and rg, obtained using various approaches. 





IHT 


HT 


d = 3 exp. 


e-exp. 


CRG 




19.13(10) 


19.31(14), 19.27(11) II 


19.06(5) 


19.55(12) [|l,pl 


22.35 |55.64| 






19.34(16) II 


19.06 11 






re 


1.86(4) 


2.1(6) H 


1.880 g 


1.867(9) (81^851 


1.74 [S6| 








1.884(32) H 








0.6(2) 




0.975 @ 


1.0(6) 1311 


0.84 II 



In order to estimate the universal quantities and r2j from the corresponding improved 
HT expansions (see App. p3 1|) , we essentially used the analysis described in Ref. Here, 
we report only the final estimates: 

(74 = 19.13(8)[2], (40) 
r6 = 1.86(3)[l], (41) 
r8 = 0.60(15) [5], (42) 

where the error in parentheses is related to the spread of the approximants and the second 
one in brackets to the uncertainty on A*, evaluated as before. Moreover, we obtained the 
rough estimate rio = —15(10). In Table ^ we compare our results (denoted by IHT) with 
the estimates obtained using other approaches, such as HT expansions for the standard 
0(3)-vector model (HT), field-theoretical fixed-dimension perturbative expansions {d = 3 
exp.), e expansions (e-exp.), and approximate solutions of continuous renormalization-group 
equations (CRG). All estimates are in good agreement. Only the e-expansion estimate of 



is significantly higher than the IHT estimate (as already noted in Ref. ||8T[, the error may 
be underestimated). The CRG estimates are much less precise than the results of other 
methods. 



B. Parametric representations of the equation of state 

In order to obtain approximations of the equation of state valid in the whole critical 
region, we use parametric representations that implement the expected scaling and analytic 
properties. We write [p7|-|89[| 



M = moR^m{e), 

t = R{l-e'^), 
H = hoR^^h{e), (43) 

where ho and mo are normalization constants. The variable R is nonnegative and measures 
the distance from the critical point in the (t, H) plane, while the variable 6 parametrizes 
the displacement along the lines of constant R. The functions m{6) and h{6) are odd and 
normalized so that m{6) = 6 + 0{9^) and h{9) = 9 + 0{9^). The smallest positive zero 
of h{9), which should satisfy 6*0 > 1, corresponds to the coexistence curve, i.e., to T < Tc 
and H ^ 0. The parametric representation satisfies the requirements of regularity of the 
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equation of state. Singularities can appear only at the coexistence curve (due, for example, 
to the logarithms discussed in Ref. [^), i.e., for 9 = 9q. The mapping (|43|) is not invertible 
when its Jacobian vanishes, which occurs when 



Y{9) = {1- 9^)m'{9) + 2p9m{9) 



0. (44) 

) are acceptable only if 9q < 9i 



Thus, parametric representations based on the mapping ([ 

where 9i is the smallest positive zero of the function Y{9). 

The functions m{9) and h{9) are related to the scaling function f{x) through 



1-9 



2 r 



X 



9^,-1 



m{9) 
■m(l) 



ni{9o) 

m{9) 
1 -s 



1 1//3 



m 

h{iy 



The asymptotic behavior (^) is reproduced simply by requiring that 

h{9) ~ {9o - 9f for 9 9q. 

The scaling function F{z) is obtained by 

z = pm{9) {l-9^y\ 
Fizi9)) = p{l-9')-''H9), 

where p may be taken as a free parameter . 



(45) 



(46) 



(47) 



C. Approximate polynomial representations 

Following Ref. |]78[, we construct approximate polynomial parametric representations 
that have the expected singular behavior at the coexistence curve (Goldstone singularity) 
and match the known terms of the small- 2; expansion of F{z), cf. Eqs. ( p6[) and (|37D. We 
consider two distinct approximation schemes. In the first one, which we denote by (A), h{9) 
is a polynomial of fifth order with a double zero at 9q, and m{9) is a polynomial of order 
(l + 2r2): 



scheme (A) : m{9) = 9 il + y^Q^^ 

V i=i 

h{9) = 9{l-9y9',)\ (48) 
In the second scheme, denoted by (B), we set 
scheme (B) : m{9) = 9, 

h{9) =9{l- 9^/9lY (1 + X^Q^^M . (49) 



i=l 
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Here h{6) is a polynomial of order 5 + 2n with a double zero at ^o- For n = the approxi- 
mations (A) and (B) coincide. Note that for scheme (B) 

Y{e) = l-e^ + 2(36^, (50) 

independently of n, so that 6i = {1 — 2(3)^^ . In both schemes, p, 6q, and the n coefficients 
Cj are determined by matching the small- 2; expansion of F{z). Thus, in order to fix the n 
coefficients q we use n + 1 values of r2j, i.e., rg, ■.■re+2n- 

As input parameters for our analysis we consider the estimates a = —0.1336(15), rj = 
0.0375(5), rg = 1.86(4), and rg = 0.6(2), which are the results of our HT analysis. The 
available estimate of rio is too imprecise for our purposes. In Fig. |I| we show the curves 
obtained in schemes (A) and (B) with n = 0, 1 and for a = —0.1336, rj = 0.0375, rg = 1.86, 
and rg = 0.6. The differences among the three approximations should give an indication 
of the uncertainty. The three approximations of F{z) are practically indistinguishable, 
and differ at most by approximately 2% (the difference between the two n = 1 curves is 
much smaller). Thus, by using the first two coefficients r2j, one obtains reasonably precise 
approximations of the scaling function F{z) for all positive values of z, i.e., for the whole 
HT phase up to t = 0. This is also numerically confirmed by the estimates of the universal 
constant F^, cf. Eq. ([38|), which is related to the large- 2; behavior of F{z). Indeed, we 
obtain F^ = 0.0262(4), 0.0266(5), 0.0266(5) respectively forn = 0, n = 1 (A), n = l (B), 
where the reported errors refer only to the uncertainty of the input parameters. This fact is 
not trivial, since the small-z expansion has a finite convergence radius Therefore, the 
determination of F{z) on the whole positive real axis from its small-z expansion requires 
an analytic continuation, which turns out to be effectively performed by the approximate 
parametric representations we have considered. 

In Fig. ^ we plot the approximations of f{x) corresponding to the schemes (A) and (B) 
for n = 0, 1, using the central values of the input parameters. The three curves are in 
substantial agreement, especially those with n = 1. Indeed, the difference between them 
is within the uncertainty due to the errors on the input parameters. These approximate 
parametric representations are not precise at the coexistence curve; indeed, as we shall see, 
the estimates of c/, cf. Eq. (|HT|), are rather imprecise and very sensitive to the value of 



rg, which is not known with high precision. We mention that in Ref. ||5^ an approximate 
expression for f{x) was obtained by approximately solving the continuous renormalization- 
group equations for the free energy (average action). The results are quite imprecise, as 
we shall show later by comparing the corresponding estimates for some universal amplitude 
ratios. 



D. Universal amplitude ratios 

From the critical equation of state one may derive estimates of several universal amplitude 
ratios. They are expressed in terms of the amplitudes of the magnetization, cf. Eqs. (pTj) 
and (PSP, of the singular part of the specific heat 

C_ff,sing = ^^1^1 (51) 

of the magnetic susceptibility in the HT phase 
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FIG. 1. The scaling function F{z) versus z. 

X = NCn-\ (52) 
of the zero-momentum four-point connected correlation function in the HT phase 



X4 



-7-2/35 



and of the second-moment correlation length in the HT phase 



(53) 



(54) 



where N — 3. We also consider the crossover (or pseudocritical) hne tmax(-f^), that is defined 
as the reduced temperature for which the longitudinal magnetic susceptibility XL{t, H) — 
dMj dH has a maximum at H fixed. The renormalization group predicts 



We consider several universal amplitude ratios: 



aA+C+ 



R 

Rc = 



52 ' 



(55) 
(56) 



(57) 
(58) 
(59) 
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fix) 




-1.0 -0.5 0.0 0.5 1.0 

X 

FIG. 2. The scaling function f{x) versus x. 



Morever, we estimate 



Ra 



R 



(C+)3' 



9a 



1/3 



T^B 
P = — 



Pr = 



Be ' 



)2(5 



R^R/^ 



R -^^ 



(60) 
(61) 
(62) 
(63) 
(64) 



Rn = 



l-Uo 



a 



(65) 



which, as suggested in Ref. |7^, should be less sensitive to the value of a than Uq. In App. 
y we give their expressions in terms of the functions m{6) and h{6). 

In Table [VI I| we report the universal amplitude ratios, as derived by the approximate 
polynomial representations of the equation of state for n = 0, 1. The reported errors are 
only due to the uncertainty of the input parameters and do not include the systematic 
error of the procedure, which may be determined by comparing the results of the various 
approximations. In Table [VI 1| we also show results for ^max, a^max and Wmax which are the 



values of the scaling variable z, x and w {w will be defined in Eq. ([7T[)) associated with the 
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TABLE VII. Results for the parameters and the universal ampHtude ratios using the scheme 



(A), cf. Eq. d 


and the scheme (B), cf. Eq. 


Note that the quantities reported in the 


first three hnes 


do not have a physical meaning, 


but are related to the particular parametric 


representation employed. Numbers marked with an 


asterisk are inputs, 


not predictions. 




n = 


n = L [A ) 


n = L [D ) 


P 


2.14(2) 


2.20(4) 


2.08(4) 




3.81(1) 


3.3(1) 


2.7(2) 


Cl 





-0.016(9) 


0.055(20) 


Uo 


1.61(2) 


1.56(3) 


1.57(3) 


Ra 


4.6(2) 


4.2(3) 


4.3(2) 




1.41(2) 


1.28(8) 


1.33(4) 


Rc 


0.173(3) 


0.19(1) 


0.184(6) 




8.2(2) 


7.7(3) 


7.9(2) 


R^ 


0.421(1) 


0.425(3) 


0.423(2) 


p 


1.201(5) 


1.17(2) 


1.18(1) 


Pc 


0.354(4) 


0.357(5) 


0.357(5) 


Rp 


2.026(6) 


2.020(5) 


2.021(7) 




1.278(5) 


1.275(5) 


1.275(6) 




8.9(1) 


8.3(4) 


8.5(2) 




1.64(2) 


1.53(6) 


1.57(4) 


7700 


0.0262(4) 


0.0266(5) 


0.0266(5) 


rs 


0.23(5) 


*0.6(2) 


*0.6(2) 


rio 
/? 


-1.1(3) 
1.28(1) 


-6(2) 
1.36(5) 


-7(3) 
1.33(3) 




0.181(5) 


0.21(2) 


0.20(1) 




-0.102(1) 


-0.105(2) 


-0.094(5) 




22(5) 


5(3) 


8(3) 



crossover line, F^, cf. Eq. (|38|) , which is related to the large-z behavior of F{z), rg and 
rio, cf. Eqs. ( p6D and (pTf), which are related to the small-z expansion of F{z), f^, and 
/g, cf. Eq. (|2^), which are related to the expansion at x = of /(x), and c/, cf. Eq. (pID, 
which is related to the behavior at the coexistence curve. Note that = where fl^ is 
related to the large-x behavior of f{x), cf. Eq. (pOl). 

From the results of Table |V11| we arrive at the final estimates denoted by IHT-PR in 
Table [VII H , obtained by taking the weighted average of the results for n = 1. The error 
we quote is the sum of the uncertainty induced by the error on the input parameters and 
of one half of the difference between the two approximations with n = 1. In most cases 
these estimates include the results of the n = approximation. In Table |VIII| we compare 



our results with those obtained in other approaches and in experiments ||101|| . We mention 
that the field-theoretical estimates of Uq have been obtained from the analysis of the fixed- 
dimension expansion in the framework of the minimal renormalization without e-expansion 
pT| , |93[| , and from the standard e-expansion to O(e^) [0. The estimate of Uq by CRG was 



obtained using Eq. (|3^) and the approximate expression for f{x) reported in Refs. |l55| , |64 
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TABLE VIII. Estimates of universal amplitude ratios obtained using different approaches. 
The numbers marked by an asterisk have been obtained by us using the results reported in the 
corresponding references. 





IHT PR 


d = 3 exp. 


e- exp. 


CRG 


HT 


experiments 




1.56(4) 


1.51(4) |91[ 

1 ^ A A ririQll 

1.544 [pJ^ 


1.521(22) [|92[ 


1.823 55,64| 




1.50(5) |7| 
1.27(9) ^ 
1.4(4) 




4.3(3) 


*4.4(4) H 
*4.46 H 


4.56(9) @ 


*3.41 pig] 

IJ PI -11 








1.31(7) 




1.33 PI 


1.11 [|5|,|4| 






Rc 


0.185(10) 


0.189(9) (9| 

0.194 g 


0.17 13 








R4 


7.8(3) 












Rt 


0.424(3) 


0.4347(20) 11 
0.4319(17) |100 


0.42 |99| 




0.431(5) m 
0.433(5) m 





See, e.g., Ref. ||102|| for a more complete review of theoretical and experimental estimates of 
universal amplitude ratios. 

In addition, from the approximate parametric representations of the equation of state, 
we obtain the estimates 

fi = 1.34(5), 
/° = 0.20(2), 
/3° = -0.10(1), 
0.0266(5), 
-6(3). 



TTICO 

^0 



(66) 



The estimate of rio should be compared with the much less precise HT result rio = —15(10) 



obtained in Sec. [IV A| . Concerning the quantities involving amplitudes at the crossover line, 
we report the estimates 



P„ = 1.18(2), 
= 0.357(5), 
Rp = 2.020(6), 
z„,ax = 1.275(5). 



(67) 
(68) 
(69) 
(70) 



In order to determine the behavior of the longitudinal magnetic susceptibility xl = dM/dH 
as a function of t and H, one may consider the scaling function 



D(w) = B-^H^-^'^ 



XL 



fix 



.1-1/5 



w 



6f{x)-^xf'ix) 



(71) 
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u 

FIG. 3. The scaling function C{u), cf. Eqs. (|7|) and @. 



The function D{w) has a maximum for t^max = 1.55(6). In order to simplify possible 
comparisons, it is convenient to consider the rescaled function 



-D(Wmax)' 

n = (72) 



which is such that the maximum corresponds to m = 1 and satisfies C(l) = 1. In Fig. |^ 
we plot the scaling function C{u) versus m, as obtained from the n = 0, 1 approximate 
parametric representations. 



E. Comparison with the experiments 

In spite of the large number of experiments, at present there is no accurate quantitative 
study of the equation of state in the critical regime. We shall discuss here three different 
representations that are widely used in the experimental work and we shall give explicit 
formulae for them. 

A first possibility ||103|| consists in studying the behavior of h/m = H\t\~'^/M versus 



= M^|t| Such a function can be easily obtained from our approximations for f{x), 
since = and 

— = k\xnf(x) (73) 
m 

where the constant k can be written as 
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k={B, 



CJ ^ 



(74) 



The universal ratio has been computed in the previous section, R^ = 1.31(7), and B 
and Be are nonuniversal amphtudes defined in Eqs. (|28|) and (p7|) . A plot of m? / B"^ versus 
C^h/m is reported in Fig. ^. It agrees qualitatively with the analogous experimental ones 
reported, e.g., in Refs. |lT7| , pO| , p!0^ . Often, for small h/m one approximates the equation of 
state by writing 



h , 2 

— = a± + o±m 
m 



(75) 



where a± and b± are numerical coefficients depending on the phase. Such an approximation 
has a very limited range of validity. In the HT phase, we obtain for —>■ 



h 
m 



Ra m? 



Rl r2n+2 f rn^ 
(2n + 1)! \B^ 



6 52 

n=2 

1 + 1.30(5)^ + 0.94(8) f^V + 0.06(2) f + 



52 



52 J 



V 52 y 



(76) 



where we have used the estimate of R4 reported in Table Vlll and the estimates of rg and 
rg reported in Table From Eq. (^) we see that the approximation (|7^) is valid only for 



very small m'^, i.e. at the 1% level only for m"^ < 0.01-02. r^j^g quadratic approximation- 



i.e. the approximation with an additional (m2)2 term- 
because of the smallness of the coefficient of m^. 



-has a much wider range of validity 
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In the low-temperature phase, Eq. (^) is theoretically incorrect, since for w? jB"^ 
we have 

h kcf / m^\^ 



1-1^ ' (77) 



m 4/?2 \^ 52 y 

where c/ is defined in Eq. (|3T1) and can be estimated roughly from the results reported in 
Table |Vli| . Eq. ( [77D is inconsistent with the approximation ( |75D near the coexistence curve, 
due to the presence of Goldstone modes. It would be correct only in Ising systems. 
Finally, note that for large we have 

A second form that is widely used to analyze the experimental data is the Arrott-Noakes 



^ ) =at + bM'/f^, (79) 



105|| scaling equation 

1/7 

where a and b are numerical constants. This approximation is good in a neighborhood of 
the critical isotherm t = 0. Since 

(80) 

using Eq. ( pQ]) and the numerical values reported in Eq. (p6|), we obtain 

J fc-V^ = (^-J + 0.96(4) t - 0.04(2) j - 0.02(2) i^- 



(81) 



Thus, at a 1% level of precision the Arrott-Noakes formula is valid approximately for 
t{MB^^)^^/^ < 25 which is quite a large interval. 

Finally, Ref. [§] reports an experimental study of the behavior of the critical system 
at the crossover line, and shows a plot of the curve C(m), cf. Eq. ([7^), in terms of the 
unnormalized variable Wexp = tH~^^'^^^\ We can attempt a quantitative comparison with 
the results reported in their Fig. 4. For this purpose, in Fig. |^ we plot C{u) in terms of 
their variable for the range of Wexp accessible to the experiment ||106|| . A direct comparison 



of this figure with Fig. 4 of Ref. H shows a very nice quantitative agreement. 



APPENDIX A: MONTE CARLO SIMULATIONS 

In this appendix we present some details of the analyses of the MC data. Details on the 
simulation can be found in Ref. W^- 



21 



0.0. 



-4.5 



-4.0 



L 

-3.5 

Log u 



-3.0 



-2.5 



exp 



FIG. 5. The scaling function C{u) versus the experimental scaling variable Uexp = tH ^/(^^). 



1. Definitions 



In all our work, considerable importance is played by dimensionless ratios (or phenomeno- 
logical couplings) R. In order to have better control on systematic errors we have studied 
four dimensionless ratios. We first consider the Binder cumulant t/4 | |107| ] and its general- 
ization Uq defined by 



where 



2i 



m 



V ^ 



(Al) 



(A2) 



is the magnetization of the system. The third quantity that we studied is the second-moment 
correlation length divided by the linear extension of the lattice ^2nA/ L. The second-moment 
correlation length is defined by 



^2nd 



XIF-I 
4 sin(7r/L)2 



(A3) 



where 



X 



He 



(A4) 
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is the magnetic susceptibility and 



^exp 



.2ttxi\ 



(A5) 



is the Fourier transform of the two-point correlation function at the lowest nonvanishing 
momentum. In order to reduce the statistical error, we averaged the results of all three 
directions of the lattice. 

The fourth quantity is the ratio Za/Zp, where Za is the partition function with anti- 
periodic boundary conditions in one of the three directions and Zp the corresponding one 
with periodic boundary conditions in all directions. Anti-periodic boundary conditions mean 
that the term J2<xy>'t'x " (py in the Hamiltonian is multiplied by —1 for x = {Li,X2,X3) 
and y = (1, 0:2, 0:3). This ratio can be measured with the help of a variant of the cluster 



algorithm, the boundary-flip algorithm. It was introduced in Ref. ||108| j for the Ising model 
and generalized to 0(A^)-invariant nonlinear a models in Ref. ||109|| . As in Refs. [^S],^, we 



use a version of the algorithm that only measures Za/Zp and does not perform the flip to 
anti-periodic boundary conditions. For a recent discussion of the algorithm, see Ref. ||36|| . 



2. Determination of R* 

First, we compute (3c and the fixed-point value of the dimensionless ratios R* for A = 
4.5 — our best approximation to A* — using the standard cumulant crossing method of Binder. 

For A = 4.5 we solve Eq. (0), computing R{L, [3) by using its Taylor expansion up to the 
third order: 

(A6) 

Here I3s is the value of (3 at which the simulation was performed, and i?, di, d2, and d^ are 
determined in the MC simulation. 



As an example, results for Za/Zp are given in Table In the fits, we include all data 
with Lmin < L < Lmax- For Ljnax = 96, x^/d-o.f. (d.o.f. is the number of degrees of freedom 
of the fit) is smaller than 1 starting from Lmin = 16. Moreover, the result for l3c is stable 
when further data are discarded. To be on the safe side, we take our final estimate from the 
fit with Lmin = 28 and Lmax = 96. 

The systematic error due to corrections to scaling is estimated by comparing the results 
corresponding to Lmin = 28 and L^ax = 96 with those with Lmin = 14 and Lmax = 48. 
We suppose that the error of Pc is proportional to L^^/'^^^ ^ L^^'^, where we assume 
pessimistically leading and not subleading corrections. Hence, we estimate the error on 
our final result as the difference between the Lmin = 28, L^ax = 96 result and the Lmin = 
14, Lmax = 48 result multiplied by 2~'^-'^/{l - 2"^-^). The systematic error of {Za/Zp)* is 
estimated in a similar fashion, assuming that the error is proportional to . 

In the same way we analyze our data for the remaining three dimensionless ratios. Our 
results are reported in Table Note that the four results for /3c are consistent within error 



23 



TABLE IX. 


Fits with ansatz 


IM\ nf niir riata 

KJl. \J U.1 U.CX L CX 


for Za/Zp at A = 4.5. We included all data with 


Lmin < L < Lmax in the fit. 












X /d.o.f. 


Pc 


[Za/Zp) 


12 


96 


5.75 


0.6862428(6) 


0.19408(3) 


14 


96 


1.90 


0.6862413(7) 


0.19419(3) 


16 


96 


1.29 


0.6862406(7) 


0.19424(4) 


18 


96 


0.62 


0.6862400(7) 


0.19430(4) 


20 


96 


0.68 


0.6862400(8) 


0.19430(5) 


22 


96 


0.75 


0.6862400(8) 


0.19430(6) 


24 


96 


0.81 


0.6862397(10) 


0.19434(8) 


28 


96 


0.60 


0.6862390(10) 


0.19443(10) 


12 


48 


5.96 


0.6862459(10) 


0.19399(4) 


14 


48 


1.92 


0.6862432(10) 


0.19413(4) 



bars. The statistical error of jSc obtained from Za/Zp and ^2nd/L is considerably smaller than 
that from U4 and f/g. As our final estimate we take /3c = 0.6862385(20), which is consistent 
with all four results. 

3. Determination of A* 

In order to compute A* we solve the equation c(A) = 0, where c(A) is defined in Eq. (^. 
In practice, we replace c(A) with its first-order Taylor expansion around A = 4.5 and thus 
evaluate A* from 

) ' • (AT) 

A=4.5/ 

In order to compute c(4.5) we fit our data for R with the ansatz (^, where we fix u; = 0.8. 
We have checked that the final result for A* has a very weak dependence on the value of a; 
used in the analysis. If we choose u) = 0.75, the results vary much less than the quoted error 
bar, indicating that the error on u) can be neglected. 

As an example, the results for f/4 at {Za/Zp) f = 0.1944 are given in Table ^ We 
see that there is a slight drift towards larger values of c(4.5) as Lmin is increased. The 
final result corresponds to Lmin = 16 and Lmax = 96. Systematic errors due to subleading 
corrections to scaling are estimated by comparing with the results obtained for Lmin = 8 
and Lmax = 48. Assuming the errors on c(4.5) to decrease as /.-'^a+i^ ~ L~^'^, we arrive at 
c(4.5) = 0.0010(6) [19], where the systematic error is quoted in brackets. 

In Table ^ we give the analogous results for U4 at {C2nd/L)f = 0.5644. Here, we see a 
larger change of c(4.5) when Lmin is varied. Also, x^/d-o.f. is larger for Lmin < 14. Since 
corrections are larger than above, we take the final estimate from the fit with Lmin = 24 
and Lmax = 96. We arrive at the estimate c(4.5) = 0.002(2) [5]. In a similar way we arrive 
at c(4.5) = 0.007(5) [15] for Uq at (Za/Zp) f = 0.1944 and c(4.5) = 0.003(4) [18] for Uq at 
(6nd/L)/ = 0.5644. 



A* ^ 4.5 - c(4.5) 
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TART F V 




dilQ C( '^i.O ) ilUlIl Lilc liL 


(|D|j Ol (74 at yZja/^pjf 




A — 4t.O. 












^max 


X /d.o.f. 


= 

R* 


c(4.5) 


6 


96 


7.87 


1.13931(2) 


0.00073(12) 


8 


96 


1.71 


1.13944(2) 


-0.00036(17) 


10 


96 


1.82 


1.13942(3) 


-0.00022(25) 


12 


96 


1.65 


1.13937(4) 


0.00039(36) 


14 


96 


1.25 


1.13930(5) 


0.00118(47) 


16 


96 


1.33 


1.13932(5) 


0.00097(60) 


18 


96 


1.33 


1.13928(6) 


0.00145(73) 


20 


96 


1.27 


1.13923(7) 


0.00217(89) 


24 


96 


1.54 


1.13920(10) 


0.00273(155) 


8 


32 


1.00 


1.13947(3) 


-0.00060(19) 


8 


48 


1.26 


1.13945(2) 


-0.00045(18) 


12 


48 


1.20 


1.13939(4) 


0.00018(38) 


16 


48 


0.75 


1.13936(6) 


0.00056(67) 



TABLE XI. 


Estimates of R* 


and c(4.5) from the fit 


(§ of ?74 at {i2nd/L)f 


= 0.5644. Here 


A = 4.5. 












^max 


xVd.o.f. 


R* 


c(4.5) 


8 


96 


9.98 


1.13984(2) 


-0.00493(19) 


10 


96 


5.05 


1.13966(3) 


-0.00313(31) 


12 


96 


2.32 


1.13948(4) 


-0.00120(39) 


14 


96 


1.30 


1.13937(5) 


0.00010(51) 


16 


96 


1.40 


1.13936(6) 


0.00022(65) 


18 


96 


1.21 


1.13929(7) 


0.00112(81) 


20 


96 


1.24 


1.13925(8) 


0.00167(99) 


24 


96 


1.50 


1.13923(11) 


0.00209(173) 


8 


48 


10.41 


1.13987(3) 


-0.00514(19) 


12 


48 


1.99 


1.13952(4) 


-0.00153(42) 



Next, we compute dc/dA at A 
finite differences 



4.5. To estimate the derivative of c(A), we consider the 



dc 
dA 



A=4.5 



[i?(L,5.0) -i?(L,4.0)] 



(A8) 



The results for our four choices of R are given in Table pCII| . We see that the results, as 
functions of L, are constant within error bars. This nicely confirms the exponent uj = 0.8. 
The final result for the derivative is obtained by averaging the results for L > 12, see 



Table PCIII] . In order to estimate the discretization error, we additionally compute the 
derivative, using R{L, A) for the pair A = 5.0 and A = 4.5. The difference with the above- 
reported result is small and in practice negligible, approximately 12% for at {Za/Zp)f = 
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TABLE XII. 


Estimates of \R{L, 5.0) 


- i?(L, 4.0)1 L'^ with CO 


= 0.8. 


L 


Ua at {Za/Zj^f 


U4 at {^2nd/ L)f 


Ue at {Za/Zp)f 


Ue at (6nd/^)/ 


6 


-0.00957(17) 


-0.01128(18) 


-0.03141(53) 


-0.03674(54) 


8 


-0.00951(21) 


-0.01136(22) 


-0.03111(67) 


-0.03697(70) 


10 


-0.00973(25) 


-0.01174(27) 


-0.03164(80) 


-0.03799(85) 


12 


-0.00972(28) 


-0.01160(30) 


-0.03150(88) 


-0.03747(93) 


14 


-0.00983(33) 


-0.01188(36) 


-0.03201(105) 


-0.03855(113) 


16 


-0.00972(37) 


-0.01181(40) 


-0.03183(115) 


-0.03848(123) 


18 


-0.00972(47) 


-0.01195(53) 


-0.03177(148) 


-0.03884(162) 


20 


-0.01022(57) 


-0.01174(63) 


-0.03392(179) 


-0.03880(194) 


22 


-0.00947(68) 


-0.01176(74) 


-0.03074(213) 


-0.03795(231) 


24 


-0.00962(79) 


-0.01209(86) 


-0.03128(249) 


-0.03915(272) 


28 


-0.00962(124) 


-0.01159(135) 


-0.03208(388) 


-0.03832(423) 



TABLE XIII. Results for [R{L, 5.0) - R{L, A)] /(5.0 - A) L'^ with lo = 0.8 for L > 12. In the 



first row A = 4.0, in the second one A = 4.5. 



C/4 at {Za/Zp)f 


U4 at {^2nd/L) f 


Ue at {Za/Zp)f 


Uq at {^2nd/ L)f 


-0.00976(16) 


-0.01177(17) 


-0.03184(49) 


-0.03823(53) 


-0.00872(27) 


-0.01035(29) 


-0.02838(86) 


-0.03356(92) 



0.1944 and f/4 at (6nd/i^)/ = 0.5644, and approximately 14% for Uq at (^a/^p)/ = 0.1944 
and Ue at (6nd/^)/ = 0.5644. 

Inserting our numerical results for c(4.5) and dc/dX into Eq. (|A7|) we get A* = 4.6(4), 
4.7(8), 4.7(8), and 4.6(8) from U4 at {Za/Zp)f = 0.1944, t/4 at {^2nd/L)f = 0.5644, Uq at 
{Za/Zp)f = 0.1944 and f/g at (^nd/-^)/ = 0.5644, respectively. The errors take into account 

the uncertainty of c(4.5), 4^ , and uj. 

aA x=A.5 

As our final result we quote A* = 4.6(4) from U4 at {Za/Zp)f = 0.1944. 



4. Critical exponents 



We compute the critical exponents u and rj using standard FSS methods. 



a. The exponent v 



The exponent v is determined by fitting the data with Eq. (^. We study the derivative 
of all four quantities [/4, U^^ Snd/-^', and Za/Zp, and fix [3f by using either (^2nd/-^)/ = 0.5644 
or (Z,/Zp)^ = 0.1944. 



As typical examples, we give fit results for 



d{Za/Zp) 



and 



dp 



/3/ 



in Tables 



PC1V| , pCV| , and pCVl| , respectively. In all these three cases, we have fixed (3f by {Za/Zp)f. 
Fixing Pf by (^2nd/-^)/ leads to similar results. 



26 



TABLE XIV. 
(ZjZp)j = 0.1944. 



Estimates of v from the fit of 



/3/ 



with ansatz (^). /3j is fixed by 



xVd.o.f. 



6 
8 

10 
12 
16 
20 



96 
96 
96 
96 
96 
96 



1.36 
0.72 
0.70 
0.71 
0.79 
0.86 



0.71215(16) 
0.71174(20) 
0.71179(26) 
0.71160(34) 
0.71137(49) 
0.71167(66) 



TABLE XV. Estimates of v 
(ZjZp)/ = 0.1944. 



from the fit of 



/3/ 



with the ansatz /3/ is fixed by 







xVd.o.f. 


V 


6 


96 


69.37 


0.70660(6) 


8 


96 


11.95 


0.70837(8) 


10 


96 


3.28 


0.70918(11) 


12 


96 


1.91 


0.70969(15) 


16 


96 


1.40 


0.71009(22) 


20 


96 


1.21 


0.71054(30) 


24 


96 


1.46 


0.71044(46) 


28 


96 


1.59 


0.71071(58) 



In the case of ^ 



/3/ 



we see that x^/d.o.f. is close to one, even if all lattice sizes L > 6 



are included in the fit. Also, the result for v stays rather stable when Lmin is increased. 



On the other hand, for 



and 



dp 



/3/ 



X^/d.o.f. comes close to one only 



starting from Lmin > 16. Also, the results for v strongly change, when Lmin is increased. It 



is interesting to notice that the estimate of v is decreasing for 



dp 



for^(%^ 



while it is increasing 



Pf 



when Lmin is increased. Assuming that this is already the asymptotic behavior, 



we can find lower and upper bounds for v. 

Taking into account the fit results for L^am > 22 we arrive at the final estimate v = 
0.7113(10). Here, the error bar includes both the statistical and the systematic error. 

Finally, we try to determine the effect of leading corrections to scaling on our estimate 
of z/. For this purpose we fit our data up to Lmax = 28 at A = 4.0, A = 4.5, and A = 5.0 
with the ansatz (^. In Table PCVII] we give our results for the derivative of the Binder 
cumulant with respect to P at {Za/Zp)f = 0.1944. In particular, for small Lmin, "we see a 



clear dependence of the result for u on A. For instance, for 
the result for A = 4.0 and A = 5.0 is 0.00176(46). 



8 the difference between 



In Table PCVIII| we give the corresponding analysis for Za/Zp at {Za/Zp)f = 0.1944. In 
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TABLE XVI. 


Estimates of 


u from the fit of 


with ansatz (S) 


/3j is fixed by 


{Za/Zp)f 


= 0.1944. 
















xVd.o.f. 




6 
8 

10 
12 
16 
20 
24 
28 




96 
96 
96 
96 
96 
96 
96 
96 




68.10 
18.76 
7.72 
4.35 
1.32 
1.43 
1.27 
1.39 


0.71826(8) 

0.71613(11) 

0.71489(14) 

0.71394(18) 

0.71261(27) 

0.71246(37) 

0.71169(56) 

0.71196(70) 


TABLE XVIL Estimates of 
where j3f is fixed by {Za/Zp)f = 


V computed from the derivative of the Binder cumulant at (3f, 
0.1944. 






-^Tnax 




xVd.o.f. 


V 


A = 4.0 


6 
8 
10 
12 




28 
28 
28 
28 




0.91 
0.75 
0.83 
1.07 


0.71301(24) 
0.71269(33) 
0.71296(44) 
0.71275(64) 


A = 4.5 


6 
8 
10 
12 




28 
28 
28 
28 




1.64 
0.74 
0.67 
0.74 


0.71226(17) 
0.71181(23) 
0.71193(32) 
0.71170(46) 


A = 5.0 


6 
8 
10 
12 




28 
28 
28 
28 




1.27 
0.76 
0.44 
0.51 


0.71153(23) 
0.71093(32) 
0.71034(43) 
0.71014(60) 



this case we see a much smaller dependence of the results for v on A. For L^^^ = 8 the 
difference between the results for A = 4.0 and A = 5.0 is —0.00065(18). The behavior in the 
case of C,2nd/ L (which is not shown here) is much the same: the dependence of the fit result 
for 1/ on A is much smaller than for the Binder cumulant. 

Taking into account the range the lattice sizes that are used to obtain the final result 
for u, we arrive at a possible uncertainty of 0.0001 for u from the derivative of Za/Zp and of 
C,2nd/L due to the uncertainty in A*. The systematic error is clearly dominated by subleading 
corrections. Our final MC estimate of u is 

z/ = 0.7113(11). (A9) 
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TABLE XVIII. Estimates of v computed from the derivative of Za/Zp at where (3f is fixed 
by (ZjZp)/ = 0.1944. 









xVd.o.f. 




A = 4.0 


8 


28 




8.97 


0.70766(13) 


10 


28 




1.63 


0.70871(18) 


12 


28 




1.58 


0.70902(27) 


A = 4.5 


8 


28 




14.38 


0.70804(9) 


10 


28 




2.84 


0.70890(13) 


12 


28 




1.57 


0.70941(19) 


A = 5.0 


8 


28 




8.52 


0.70831(13) 


10 


28 




2.89 


0.70916(18) 


12 


28 




2.40 


0.70965(25) 



TABLE XIX. Estimates of rj from fits of the magnetic susceptibility at A = 4.5 with Eq. (p!o|). 
is fixed by {Za/Zp) f = 0.1944. 



X /d.O.f. T] 

12 96 32.55 0.03557(5) 

16 96 6.36 0.03641(7) 

20 96 1.73 0.03682(9) 

24 96 1.00 0.03710(13) 

28 96 0.81 0.03725(16) 

32 96 0.81 0.03740(24) 



h. The exponent t] 

We determine the exponent r] by using Eq. (|l^) and also a fit ansatz that includes a 
constant background term, Eq. (|TT|). We fix Pf by setting either {C,2nd/ L)f = 0.5644 or 
(V^p)/ = 0.1944. 

Fits for A = 4.5 with the ansatz (|lO|) are given in Table XIX| (/5/ fixed by (Za/Zp) f = 



0.1944) and PCX| (/?/ fixed by {$,2nd/ L)f = 0.5644). In both cases, x^/d-o.f. becomes close 
to one starting from Lmin = 24. Moreover, in both cases the fit results for rj are strongly 
increasing as L^in is increased. For Lmin = 32 we have a consistent result of t] = 0.0374(2). 



Next, we checked the dependence of the result for r] on A. In Table PCXI] we give results 
for A = 4.0, 4.5, and 5.0 for Lmin = 16 and Lmax = 28. We see a rather strong dependence 
on A. The difference between the results for A = 4.0 and 5.0 is 0.00080(25). Taking into 
account the range of lattice sizes used to obtain our final estimate, we arrive at an error of 
0.0002 on f] due to the error on A*. 

Finally, we performed fits with the ansatz ([Tl|) . The results are summarized in Tables 



29 



TABLE XX. Estimates of i] from fits of the magnetic susceptibility at A = 4.5 with Eq. ([lO|). 
Pf is fixed by (6nd/^)/ = 0.5644. 



X^/d.o.f. T] 

12 96 33.81 0.03592(4) 

16 96 8.24 0.03656(5) 

20 96 2.35 0.03695(7) 

24 96 0.77 0.03725(10) 

28 96 0.59 0.03735(12) 

32 96 0.65 0.03742(18) 



TABLE XXL Estimates of r] from fits of the magnetic susceptibility at A = 4.0, 4.5, 5.0 with 
Eq. (loD. Pf is fixed by (Za/^p)/ = 0.1944. 

X^/d.O.f. T] 

40 16 28 475 0.03610(18) 

4.5 16 28 2.96 0.03562(13) 

5.0 16 28 8.53 0.03530(18) 



TABLE XXIL Estimates of rj from fits of the magnetic susceptibility at A = 4.5 with Eq. (|TT|). 



pf is fixed by {Za/Zp)f = 0.1944. 







xVd-o.f. 


V 


b 


8 


96 


2.18 


0.03832(8) 


-0.657(9) 


10 


96 


1.22 


0.03811(10) 


-0.617(17) 


12 


96 


0.72 


0.03790(12) 


-0.555(26) 


16 


96 


0.66 


0.03782(17) 


-0.528(54) 



XXI] and XXIIL We observe that a x /d.o.f. close to one is already reached for 



10. 



16 the results 



Moreover, the result for rj changes little with increasing Lmin- For L^i^ 
obtained by fixing Pf by Za/Zp and ^2nd/L agree. Therefore, we give as our final result 
1] = 0.0378(6). The error bar is such that it includes the result of the fits with ansatz (p^Of ). 



APPENDIX B: ANALYSIS OF THE HIGH-TEMPERATURE EXPANSIONS 



In this appendix we report a discussion of our HT analyses. It should allow the reader 
to understand how we determined our estimates and the reliability of the errors we report, 
which are to some extent subjective. More details on the methods we use are reported in 
Ref. M 



30 



TABLE XXIII. Estimates of rj from fits of the magnetic susceptibility at A = 4.5 with Eq. (11). 
/?/ is fixed by (6nd/^)/ = 0.5644. 



X^/d.o.f. 7] b 

8 96 (180 0.03756(6) -0.379(6) 

10 96 0.68 0.03765(7) -0.396(11) 

12 96 0.64 0.03772(9) -0.416(20) 

16 96 0.53 0.03780(13) -0.454(39) 



1. Definitions and HT series 



We computed the HT expansion of several quantities for the 0^ lattice Hamiltonian (|l|) for 
generic values of A by using the linked-cluster expansion technique. A general introduction to 
this technique can be found in Refs. ||110| - |112(] . We calculated the 20th-order HT expansion 



of the magnetic susceptibility and of the second moment of the two-point function, 

X = rn^ = 5^x2(0„(O)0„(x)), (Bl) 

X X 

and therefore, the second- moment correlation length = 7712/ (6x)- Moreover, we computed 
the HT expansion of the zero-momentum connected 2j-point Green's functions X2j 

X2j = 5Z (<^"i(O)<^"i(^2)...0a,(a;2i-i)0a,(a;2i))c (B2) 

X2,---,X2j 

{x = X2)- More precisely, we computed Xi to 18th order, xe, Xs, and Xio to 15th order. In 
Table pCXlV| we report the series for the 0^ Hamiltonian with A = 4.5. We chose this value 



because it is very close to the best estimate of A*, and because for this value of A we have a 
precise MC estimate of pc, Pc = 0.6862385(20). 

The HT series of the zero-momentum four-point coupling and of the coefficients 
that parametrize the small-magnetization expansion of the equation of state can be computed 
using their definitions in terms of X2j and i.e.. 



and 



re = 10 -4^^, (B4) 
3(iV + 4) xl ^ ^ 



= 2Sn - 280(iV + 2) X6X2 35(iV + 2)^ xsxj 
3(iV + 4) xl 9(iV + 4)(iV + 6) xl ' 
r = 11^^(11+^X6X2 mN + 2Yxlxl 



+ 



(iV + 4) xl (iV + 4)2 xt 
1400(iV + 2)2 xsxl 35(iV + 2)3 xioxl 



3{N + 4){N + 6) xl S{N + A){N + Q){N + 8) xt ' 
The formulae relevant for the Heisenberg universality class are obtained setting = 3. 
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TABLE XXIV. CoefHcients of the HT expansion of m2, x, X4, Xejj XSi and xio, for the cf)^ 
Hamiltonian with A = 4.5. 



m2 



X2 



Xi 





1 0.61164859624923922306340 

2 2.34346567036967995931352 

3 6.12631563804734885224064 

4 13.7411678013148223984708 

5 28.1806907736271607364126 

6 54.6110112231156495979858 

7 101.601106321810555490226 

8 183.444896900967255568497 

9 323.515763329983708761073 

10 560.008592676223985571196 

11 954.596109677157386424652 

12 1606.62254117132464855356 

13 2674.82852285795544112124 

14 4412.16935517377929559254 

15 7219.36840082433737483629 

16 11729.5598928466629681760 

17 18938.5632727811981324466 

18 30408.8505644977121915813 

19 48583.1353446096892196268 

20 77271.9486733817666941548 



0.95784805390722532625540 
1.83494578874771766919200 
3.02570526555699314952643 
4.91084272357011167693574 
7.72413685583625911405950 
12.0696651111374510782698 
18.5898743639602203433897 
28.5129469684875374524852 
43.3910414400431056256283 
65.8457097219033943862880 
99.4318415601856883686915 
149.842255361185424824490 
225.053071205898843291857 
337.491827618915172447234 
504.872666999718714410906 
754.353044416651842398902 
1125.02916338959766441444 
1676.21770934179961650320 
2493.83291987123696064787 
3707.29637073719901187394 
5504.79157669035824056791 



-0. 10220686631889066464185 
-0.7831891839960428631 1017 
-3.40993303251195614433724 
- 1 1 .4547261 149387426 1 1 1654 
-32.8622857609807070588938 
-84.8703770522664536068131 
-203.048391944861326183480 
-458.487070991901596469076 
-989.113871576954399902821 
-2056.80132192600410402482 
-4148.88336851312159356144 
-8158.09343927729005801740 
-15696.1032586702224741180 
-29637.6862646277928041257 
-55053.5845813513443015904 
- 100803. 117025457078176935 
-182227.804691700019121686 
-325689. 103869085626575337 
-576156.574987690512078391 



X6 



X8 



Xio 



0.168561977829196181908019 

1 2.564255760481036314377282 

2 19.75833949883362663601629 

3 107.6153442827461300177955 

4 470.3177686849044939915090 

5 1762.676494500452240334969 

6 5892.468315770018895075220 

7 18026.10606034742588181350 

8 51364.40898933529517266700 

9 138079.9439575146969655288 

10 353553.4243627006337358070 

11 868613.6815956677724696610 

12 2059460.213048375196521201 

13 4734189.661963741454010320 

14 10591072.53372872753915110 

15 23130642.47447337362164738 



-0.61605090918722894152663 
-15.2300029871348773005564 
-178.372137615950383513749 
- 1404.03404775262637194221 
-8533.18422765821528578124 
-43117.7286797414375624292 
- 189462.996356270457135473 
-745826.772293316301527782 
-2685696.43890603392620130 
-8982884.52144712077544905 
-28231531.3374765513159336 
-84122578.3791382746491922 
-239356596.354404154637789 
-654084831.403386579980857 
- 1724768015.94141791734955 
-4405965912.50799258526090 



3.90323863425838698945752 
141.589060417753517915167 
2318.78554766886373709723 
24643.8854034121990204490 
196815.472836825966956507 
1278183.93414715336465143 
7085909.43897471702725255 
34639336.4985891651938076 
152792638.179196536855095 
618543549.844101982857471 
2328085189.97290506149579 
8229819042.53291083435060 
27546203927.2418206714112 
87876514492.4881461446469 
268647512699.595417711632 
790611696518.867100914318 
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2. Critical exponents 



In order to estimate 7 and we analyzed the 20th-order HT expansion of the magnetic 
susceptibihty and the 19th-order HT expansion of f5. We analyzed the HT series by means 
of integral approximants [|113|| (lA's) of first, second, and third order (lAl's, IA2's and IA3's 
respectively). Since the most precise results are obtained by using the MC estimates of Pc 
to bias the approximants, we shall only report the results of the biased analyses. We used 
the values of Pc obtained in Sec. [A 2|, i.e. 



p^{X = 4.0) = 0.6843895(35), (B5) 
/3c(A = 4.5) = 0.6862385(20), (B6) 
p^(X = 5.0) = 0.6875638(37). (B7) 

We considered several sets of biased lA's, and for each of them we obtained estimates of the 
critical exponents. In the analysis we followed closely Ref. I^Uj. Thus, in the following we 
shall heavily refer to it for notations and a more detailed description of the analyses. 

Given an nth-order series /(/5) = X^iLo^i/?*, its fcth-order integral approximant 
[mfc/mfc_i/.../mo//] lAA; is a solution of the inhomogeneous fcth-order linear differential 
equation 

Pkmf^'\f3) + p,.mf'-'\(^) + - + A(/5)/(^)(/3) + pmf{f^) + R{f^) = 0, (B8) 

where the functions Pi{P) and R{f3) are polynomials of order rrii and / respectively, which 
are determined by the known nth-order small-/3 expansion of /(/?). We considered two types 
of biased lA/c's: 

(i) The first type of biased lA/c's, which will be denoted by blA/c's, is obtained by setting 

Pk{P) = il-P/Pc)PkiP), (B9) 

where Pk{(3) is a polynomial of order — 1. 

(ii) Since on bipartite lattices (3 = —(3c is also a singular point associated to the antifer- 



romagnetic critical behavior ||114|| , we consider lA/c's with 



PM={^-^'/PI)pM, (BIO) 

where Pk{f3) is a polynomial of order — 2. We shall denote them by h±lAk's. 

In our analyses we considered diagonal or quasi- diagonal approximants, since they are 
expected to give the most accurate results. For each set of lA/c's we calculated the average 
of the values corresponding to all nondefective lAfc's. Approximants are considered defective 
when they have singularities close to the real jS axis near the critical point. We also discarded 
some nondefective lA's — we call them outliers — whose results are far from the average of 
the other approximants. All details can be found in the App. B of Ref. PB . 



In Tables PCXV| and pCXVIj we report the results for 7 and u respectively, obtained by 
analyzing the series for A = 4.0,4.5, and 5.0. There, we also quote the "approximant ratio" 
ra = {g — f)/t, where t is the total number of approximants in the given set, g is the number 
of nondefective approximants, and / is the number of outliers which are discarded using an 
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TABLE XXV. Results for 7 obtained from the analysis of the 20th-order HT series of x- 



A 


approximants 


ra 


7 


4.0 


blAl 


(35 - 3)/48 


1.39508(6) [32] 




bIA2 


(77- 7)/115 


1.39503(16)[33] 


4.5 


blAl 


(36 - 3)/48 


1.39585(4) [18] 




h , T A 1 


('91 1 /4R 






bIA2 


(93 - 11)/115 


1.39580(10)[18] 




b±IA2 


(84- 7)/100 


1.39579(18)[18] 




bIA3 


(56 - 6)/61 


1.39582(5) [19] 


5.0 


blAl 


(34-3)/48 


1.39652(6) [32] 




bIA2 


(107- 13)/115 


1.39648(7) [34] 


TABLE XXVL 


Results for u obtained from the analysis of the 19th-order 


series of 


A 


approximants 




V 


4.0 


blAl 


(37-6)/37 


0.71061(1)[14] 




bIA2 


(65 - 3)/70 


0.71055(17)[14] 


4.5 


blAl 


(37-4)/37 


0.71110(3) [8] 




b±IAl 


(31 - 3)/36 


71111(2^81 




bIA2 


(67- 2)/70 


0.71108(6) [7] 




b±IA2 


(54-3)/55 


0.71114(3) [8] 




bIA3 


(26 - 2)/34 


0.71110(10)[10] 


5.0 


blAl 


(36 -4)/37 


0.71151(5)[15] 




bIA2 


(67- 5)/70 


0.71154(6) [14] 



algorithm described in App. B of Ref. [0; g — f is the number of "good" approximants 
used in the analysis. Notice that g ^ f, and g — f is never too small. For each analysis, 
beside the corresponding estimate, we report two numbers. The number in parentheses, ei, 
is basically the spread of the approximants for /3c fixed at its MC estimate. It is the standard 
deviation of the results obtained from all "good" lA's divided by the square root of r^, i.e., 
Ci = o'/y/va. The number in brackets, 62, is related to the uncertainty on the value of /3c 
and it is estimated by varying /?c in the range [/?c — ^Pc, Pc + A/3c]. 

APPENDIX C: UNIVERSAL AMPLITUDE RATIOS FROM THE PARAMETRIC 

REPRESENTATION 

In the following we report the expressions of the universal amplitude ratios in terms of 
the parametric representation (|43|) of the critical equation of state. 

The singular part of the free energy per unit volume can be written as 

^sing = homoR^-''g{e), (CI) 
where g{6) is the solution of the first-order differential equation 
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(1 - e')g'ie) + 2(2 - a)egie) = Y{e)h{e) (C2) 

that is regular at 6' = 1. The function Y{6) has been defined in Eq. (44). The longitudinal 
magnetic susceptibility can be written as 

The function g2{0) must vanish at in order to reproduce the predicted behavior at the 
coexistence curve Xl ~ H"^^"^, according to 

92{0) ^00-6 for 6^ Oq. (C4) 

From Eq. ( |C3D we see that g2{0) satisfies this condition if h{9) ~ (6*0 — 9Y for 9 9q. 

From the equation of state one can derive universal amplitude ratios of zero-momentum 
quantities. We consider 

f/o = A-^/A- = {91 - If-^y (C5) 

Rx ^ = (00 - ir[m{9o)r'[m{l)]-'h{l), (C6) 

c 

Rc ^ — ^ = -«(1 - «)(2 - aM - ir^[mi9o)]-'giO), (C7) 
R.^-^ = P^mi9o)r{9l-iy'r (C8) 



Using Eqs. (^) and (|47|) one can easily derive the expressions of the various coefficients 
that characterize the asymptotic behavior of the scaling functions f{x) and F{z), such as 
c/,/° for f{x) and F°°,r2j for F{z). Concerning the ratios involving amplitudes along the 
crossover line, one finds 

Pm = ^5 = ^max/l-^max) j (C9) 
-De 

Pc = -^V^ = i^(^max)"', (CIO) 

Rp = ^ = i^'(^max). (Cll) 
Op 

Here Xmax and ^max are the values of the scaling variables x and z computed at 6'max, where 
^max is the solution of the equation 

f36F[z{9)\F"[z{9)\ - -iF'[z{9)f = 0. (C12) 
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